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Abstract 

Dark matter with mass below about a GeV is essentially unobservable in conventional direct detection 
experiments. However, newly proposed technology will allow the detection of single electron events in 
semiconductor materials with significantly lowered thresholds. This would allow detection of dark matter 
as light as an MeV in mass. Compared to other detection technologies, semiconductors allow enhanced 



r~| ■ sensitivity because of their low ionization energy around an eV. Such detectors would be particularly 

Oh' 
I ■ sensitive to dark matter with electric and magnetic dipole moments, with a reach many orders of magnitude 

Q-i' 

2h \ beyond current bounds. Observable dipole moment interactions can be generated by new particles with 



masses as great as ^ 10 TeV, providing a window to scales beyond the reach of current colliders. 



1. Introduction 
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cn 

. ' The particle nature of dark matter (DM) has been well established by astronomical and cosmological 

cn ■ 

^) ' data (||l|, g] and references therein). It is reasonable to expect the DM particle to carry non-gravitational 
interactions. These non-gravitational interactions may permit direct detection of DM, plausibly leading 
to a deeper understanding of its origins and the structure of particle physics. A good case can be made 
for weak-scale interactions between the standard model and DM. The possible existence of new states 
at the weak scale, as suggested by the hierarchy problem, could lead to such interactions. A variety of 
experiments are currently probing these interactions. These experiments measure the energy deposited by 
DM as it scatters off the atoms in a detector. Current experiments are sensitive to recoil energies > keV 
ly]. At these recoil energies, the experiments are dominantly sensitive to the scattering of DM particles 
with masses larger than ~ 1 GeV off the atomic nucleus [0]. 
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A lighter particle bound to the DM halo is kinematically forbidden from depositing energies greater 
than a keV. Owing to our ignorance of the physics responsible for DM, it is desirable to develop technological 
tools to explore all possible regions of the DM parameter space. The ability to detect low energy (~ eV) 
electron recoil events will significantly extend our reach into this parameter space, as was demonstrated 
recently in [^. The energy deposited by a light DM particle on the nucleus is suppressed by the nuclear 
mass. The difficulty of detecting such low energy nuclear recoils is further complicated by the anaemic 
response of the nucleus to such events. However, since the electron is light, DM can dump more energy into 
it. Further, energy deposition into electrons can lead to more readily identifiable events such as ionization 
in the detector. Technological advances in semiconductor-based DM detectors have made the detection of 
such ionization events a realistic possibility. Strategies to suppress backgrounds to ultra low levels similar 
to typical direct detection experiments also seem feasible ||6|, 0]. One possible experiment of this type is 
the CDMSLite proposal, which would modify existing CDMS technology to reduce the energy threshold 
by approximately three orders of magnitude, thereby allowing detection of single electron recoils |§]. 

The ionization of electrons always involves transfer of momentum from DM to the atomic nucleus. The 
cross-section for such interactions is suppressed by a form factor for momentum transfers much bigger than 
the inverse Bohr radius of the concerned electron. The electrons in semiconductors and noble gas detectors 
have comparable Bohr radii. However, the bandstructure of the semiconductor allows for the ionization of 
electrons with relatively lower energy (~ 1 eV) in comparison with the energy needed to ionize electrons 
in a noble gas detector (~ 10 eV). Since the DM particle has to lose more energy in the case of a noble gas 
detector, the momentum transferred to the nucleus is higher, leading to a form factor suppression of the 
cross-section. Owing to the smaller energy transferred in the case of the semiconductor, the momentum 
transferred to the nucleus is also smaller, leading to an unsuppressed cross-section. 

In this paper, we argue that semiconductor detectors able to measure the production of single electron- 
hole pairs have the potential to detect light DM in a wide range of parameter space, orders of magnitude 
beyond current bounds. We also show that such semiconductor devices possess an enhanced sensitivity 
to light DM in comparison with noble gas detectors. We begin, in section g, by considering the possible 
models for DM-electron interactions, either through renormalizable couplings or effective operators such as 
electromagnetic dipoles. We illustrate, through the aid of a simple concrete example, the ease with which 
these electromagnetic moments are generated for DM particles as a consequence of new states at the weak 
scale. In section H, we compute the rate for a light DM particle to scatter off a valence electron bound to 



a semiconductor. In computing this rate for various operators, we focus in particular on electromagnetic 
moments because they are both easily generated in many models and have enhanced cross-sections due to 
their coupling to a long range force carrier, namely, the photon. Consequently, these operators may offer 
the easiest way to probe the existence of such light DM. Using an estimate of the possible backgrounds 
at CDMSLite, we examine the sensitivity of such a device to DM-electron interactions in section |^, in 
comparison to current bounds on these operators, as well as the limits possible with noble gas detectors. 

2. Models 

In order to study the potential reach of detectors such as CDMSLite, we must consider the possible 
operators that generate DM-electron interactions, as well as the current constraints on the corresponding 
parameters. We restrict ourselves to the simplest extensions to the standard model, but the DM sector 
could contain other scenarios (such as |9-13|), possible signals of which should be studied in future work. 
We also focus specifically on the case of light DM (me < m^ < 10 GeV), in order to find unexplored 
parameter space for these simple models. 

2.1. Dipole Moments 

The simplest extension to the standard model is for DM and electrons to interact electromagnetically. 
Due to constraints on the possible electric charge of DM ||l^, the lowest dimensionality available for 
electromagnetic interactions corresponds to the dimension-five dipole moment operators 

Cdipoie = -^x^^^(Aix + dxl'^)xF^,u, (1) 

where ^^ and dy. correspond to the DM magnetic and electric dipole moments, respectively. These dipole 
moments correspond to a cutoff scale {d^ ~ ^^^) and arise from loop interactions involving heavy charged 
particles, such as those shown in Figure |^. These dipole moments are then easily generated in models where 
DM carries a conserved charge, such as asymmetric dark matter [15-p^]. As an illustrative example, we 



consider two new heavy intermediaries (a fermion and a scalar), in the limit where both have approximately 
the same mass M and coupling to DM g. We then obtain (based on calculations similar to ||2^) the dipole 
moment 
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Figure 1: One-loop contributions to DM dipole moment due to a charged fermion-scalar pair. 
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Figure 2: One-loop contribution to DM mass due to a charged fermion-scalar pair. 

with the same approximate form for ^u^. There might be some worry that loops involving these heavy 
charged intermediaries, shown in Figure |2|, would push the natural DM mass beyond the MeV or GeV 
scale. The contribution from this diagram is 

5m^ - ^^. (3) 

The important feature of this expression is that decreasing the coupling g between DM and the heavy 
intermediaries decreases the effective scale contributing to the DM mass. However, this decrease in g 
actually increases the effective scale contributing to the DM dipole moment. This means that for a generic 
set of heavy charged intermediaries, a large effective dipole scale does not imply a large mass contribution, 
provided the coupling with DM is small. For example, a charged fermion-scalar pair with M ~ 500 GeV 
and g ~ 0.2 would contribute 5m^ ~ 100 MeV and d^ ~ ?> x 10~^ TeV~^. As we will show in section 
0, the enhanced cross-sections of dipole interactions at low momentum transfer make them the strongest 
candidate for detection with CDMSLite, with experimental sensitivity to effective mass scales < 10'^ TeV. 

2.2. Effective Pointlike Vertex 

The next simplest extension is the dimension-six effective four-fermion vertex, which corresponds to 
the exchange of a very massive mediator (such as a scalar or vector) which is then integrated out of the 
theory. An example is the vector-channel operator 



Cpoint = ^XT^X^eTAtV'e, (4) 

where A corresponds to the cutoff of this effective theory (roughly the mass of the intermediary particle) . 
The strongest constraints on these pointlike interactions come from collider experiments such as LEP 



[21-23]. For example, the vector interaction above is currently restricted to a cutoff scale A > 480 GeV. 
Calculations based on the method presented in section y indicate that CDMSLite would only be able to 
search for pointlike interactions up to the scale A w 200 GeV. The weakness of this projected sensitivity is 
due to the lack of enhancement for DM pointlike scattering at low recoil energies. Semiconductor detectors 
will then be no more sensitive to dimension-six DM interactions than collider experiments. However, higher 
scale physics which generates dimension-six operators will very often generate dipole moment operators, 
as well, whose recoils suffer from less suppression. As will be shown in section ^, dipole operators will 
therefore allow CDMSLite to probe scales far beyond the reach possible with pointlike interactions, either 
in direct detection or collider experiments. In light of these facts, we do not consider pointlike interactions 
for the remainder of this paper. 

2.3. Broken U{1) 

The final possibility for simply extending the standard model is to introduce new light particles to 
the theory. In this case, DM can interact via a broken U{1) gauge interaction, with the corresponding 
dimension-four operator 

C-A'x = 9x^'f,Xl^X, (5) 

where A' is the DM gauge field. A DM-electron interaction could then result from kinetic mixing of 
amplitude e between A'^ and the standard model photon, which can be diagonalized to give a DM gauge 
boson-electron coupling term 

Ca'c = -eeA'^i^ej^tPe- (6) 

This interaction also has a large parameter space available for semiconductor detectors, which was 
discussed in [^. We consider this interaction later in some detail for the sake of of thoroughness and 
comparison with that work, but our main focus is on the potential exploration of new physics through DM 
dipole moments. 



3. Detection Rate 

As a simpler conceptual example, we first consider DM ionizing a single isolated atom. We then turn to 
DM interacting with a semiconducting lattice to excite an electron from the valence band to the conduction 
band. Our approaches to these two cases are rather similar, as are the resulting cross-sections. This topic 
is similar to ||5|, |2^, 25], but differs in the treatment of the lattice bandstructure and in the inclusion of 



momentum transfer to nuclei. 
3.1. Basic Kinematics 

First, we briefly review the kinematics of electron recoil interactions. Our convention for momenta is to 
use pto indicate incoming momenta and k for outgoing momenta. Also, since the DM velocity w^ ~ 10^'^, 
we use simpler nonrelativistic kinematics. 

Interactions are classified by the recoil energy Eji, defined as the kinetic energy of the outgoing electron 

In the lab frame, the electron (initially in a bound state) can have a nonzero incoming momentum pe, 
with the probability for this momentum determined by the electron's momentum space wavefunction ■(/;. 
However, the initial energy for the electron is simply —Eb, the binding energy associated with its initial 
state ip. We also choose the lab frame to be such that the nucleus has no initial momentum (pN = 0). 

The momenta values we consider are much lower (~ keV) than the nuclear masses of silicon and 
germanium, such that the final kinetic energy of the nucleus can be neglected. The resulting energy 
conservation equation can be rewritten as 

kl=pl-2m^{ER + EB). (8) 

We can then calculate the minimum possible momentum transfer q necessary to ionize an electron 
with binding energy Eb- For the case of semiconductors, with Eb ~ 1 eV, the DM must at least transfer 
momentum g ~ 1 keV for ionization to be possible. For noble gases, with a larger Eb ~ 10 eV, the mininum 
transfer necessary is q ^ 10 keV. This increase in momentum transfer away from the inverse Bohr radius 
and into the form factor regime suppresses the noble gas detection rates, reducing their experimental 
sensitivity in comparison with semiconductors. This suppression in the case of dipole moment interactions 
is given in detail in equations (|lq ) and ([iq). 
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3.2. Single Atom Ionization 

Our example initial state consists of free DM and a bound hydrogenic atom. If the atom is ionized 
by the DM-electron interaction, then the final state consists of the recoiling DM, escaping electron, and 
remaining nucleus. For calculational simplicity, we model the nuclear final state as a plane wave. However, 
the final state of the nucleus is in fact not relevant to the cross-section. Rather, what matters is that the 
nucleus, due to its large mass, can absorb momentum at negligible energy cost. From the perspective of the 
DM-electron system, the recoiling nucleus then breaks momentum conservation while preserving energy 
conservation. 

The electron's resulting wavefunction is also deformed by the presence of the charged nucleus, causing 
it to deviate from a simple plane wave. This deformation causes a substantial enhancement to the cross- 
section, which can be approximated by combining a plane wave final state with a momentum-dependent 
enhancement factor. This factor is similar to the standard treatments of beta decay [26| and Sommerfeld 



enhancement [^ (for a clear review, see p^]), as well as the work on noble gas detectors in ||5|. This 
enhancement factor can be found by exactly solving the Dirac equation for a free electron in the presence 
of a Coulomb potential, and comparing the solution to that of a plane wave, yielding 



nke) = Y^^, (9) 



where v is the fee-dependent factor 



^^Z^ffTU^^ (10) 

In order to calculate the cross-section for DM to ionize the atom, we then simply need to look at the 



usual scattering formula [29| 

""^H^fn^J^-^i^z-BOK/mr (11) 

Using this formula, we can then write out the full cross-section for our ionization process, assuming 
(as stated earlier) that the lab frame corresponds to the nucleus rest frame and using the approximation 
nie << rriN, 

F(ke) d?ky d?ke d^kiy . ,a<-4,, x t,^ r r x ^ .-^ /^ r n ^ /,^x 

d<r = -^ j^ j^ j^{27r)^6\kf - p,) i,{p^ -k^- ke) H,nt{Px - kx) , (12) 



where ip and Hint refer to the Fourier transforms of the initial electron bound state wavefunction and DM- 
electron interaction Hamiltonian, respectively. For our initial bound state, we use the hydrogen ground 
state wavefunction 



ipo{x) 
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-r/a 



(13) 



where a is the Bohr radius. 

We now calculate the cross-sections for the interactions considered earlier, the first of which is the 
electric dipole moment (EDM) interaction. Averaging over initial spins, summing over final spins, and 
using the approximation kg « q (for details, see the appendix), we obtain the approximate cross-section 
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(14) 



6(1 + aWrmnf 

The term v^ refers to the incoming velocity of the DM particle (in the lab frame), and the various 
momenta have the following definitions 



ke = yj2meER, 

Imin = m^v^ - Jm\v\ - 2m^{ER + Eb). 



(15) 



In order to understand the form factor suppression for momentum transfer above a ^, we can take the 
further limit of aqmin >> 1) obtaining 
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(16) 



The second cross-section corresponds to the similar magnetic dipole moment (MDM) interaction. Using 
the same approximations and variables as the EDM case, we find the cross-section 
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(17) 



Not surprisingly, this cross-section is quite similar to that of an EDM, but is suppressed by an additional 



approximate factor of c? . This factor corresponds to the average velocity of the bound electron ((we 
We can also obtain the similar form factor 



a). 



da 



MDM 



l6a'^filkeF{ke 



(18) 



The final cross-section corresponds to the broken f/(l) DM gauge interaction, for which we consider 
two hmiting regimes. The first corresponds to a heavy mediator, with mass m^ much greater than the 
momentum transfer, such that the interaction is effectively pointlike. In this limit, the cross-section is 



da 



heavy 



128X^keF{ke 



where A is the effective DM-electron coupling 



1 



(1 + a^qi.J^ '■ 



(19) 




(20) 

The second limit corresponds to a light mediator, with mass much less than the momentum transfer, 
such that A' is effectively massless. This yields the cross-section 
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3.3. Semiconductor Valence Band 



1 + g Qmin 



(21) 



We now consider an atom in a semiconducting lattice. For an electron in any periodic potential, the 
delocalized eigenfunctions of the Hamiltonian can be expressed in terms of localized wavefunctions, in the 
form |0|,|31| 



V'r(x) 



1 



E' 



,ib-x„ 



[X - Xr 



(22) 



where A^ is the total number of lattice sites and 6 is a wavevector with components related to the dimensions 
Lj of the lattice by the relationship 



27rni 



(23) 



where i runs over the values x, y, and z, and rij is any integer from 1 to A'j, the number of lattice sites in 
the i-direction (N = N^NyN^). The wavefunction (j){x — Xn) is a localized Wannier wavefunction centered 
around each individual lattice site (located at x„). 

In the tight-binding, linear combination of atomic orbitals (LCAO) approximation [32[, the Wannier 
wavefunctions (j) are written in the basis of free atomic orbitals. These coefficients are very small for all 
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atomic orbitals except those near the bound-state energy of i/jr |30|], which for valence band states are the 
highest occupied s- and p-states. These outermost states can be reasonably approximated with hydrogenic 
wavefunctions, due to the screening of the nuclear charge by the inner core electrons. 

For the final state, the valence band electron is excited into the conduction band, where it can be 
treated as an approximately free electron, with two corrections. Due to the weak periodic potential of 
the lattice, a conduction band electron propagates with the effective mass ml = feme- The correction 
factor fe is an element-dependent factor determined by the direction of the electron momentum in the 
lattice and the energy curvature along the conduction band. We use the density of states average values 
for /e corresponding to the very edge of the conduction band. This gives an approximate estimate for the 
interaction cross-section, which will only be slightly modified by a more exact calculation. We use fe = 1-1 



for silicon and fe = 0.6 for germanium [p3| , 34]. 

The second correction comes from the presence of the positively charged hole remaining in the valence 
band. The Coulomb interaction between these charges causes the same enhancement as the atomic case, 
where the effective charge Zejf felt by the outgoing electron is simply that of the remaining hole. 

With these particular initial and final states, we can use a similar approach to the free atom considered 
earlier, writing the 6-dependent interaction cross-section 
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(24) 



For interactions localized to a single lattice site (momentum transfer of 0{a~^)), the cross terms 
for wavefunctions of different sites are negligible. This cross-section can then be expressed in terms of 
interactions with a single local Wannier wavefunction. We therefore approximate cross-sections involving 
highly delocalized states spread throughout the entire lattice by calculating cross-sections involving a single 
localized state, repeated periodically at the A^ sites of the lattice. 
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Our full cross-section for a single lattice site is then simply the average of the individual dar, 



da = — > dar. 



(25) 



(26) 
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Figure 3; Density of states as a function of binding energy for the valence bands of germanium (red) and silicon (blue), 
normalized such that J dEs piEs) = 1- 

In this approximation, the 6-dependence for each dar is contained entirely in the initial state binding 
energy Eb- Our total cross-section then changes from an average over all possible b to an integral over all 
possible Eb, 



da ^ / dEB p{EB)da{EB), 



(27) 

where p{Eb) is an experimentally determined density of states (based on [^, ^] and shown in Figure |3|) 
accounting for the fact that some Eb values correspond to more ipr states than others. This density of 
states then serves as an efficiency factor for scattering at various binding energies. For example, there is 
zero detection efficiency in germanium at the minimum Eb of 0.7 eV, but the efficiency rapidly increases 
for slightly larger Eb- Note that our final result does not contain any directional dependence resulting 
from the lattice structure, but rather gives the directionally-averaged behavior of the total cross-section. 

The cross-sections for electrons in a semiconducting lattice are then accurately approximated by those 
of a free hydrogenic bound state, in a weighted average over initial-state binding energies, with an altered 
final-state electron mass. 

3.4- Detection Rate 

Once the interaction cross-section is known, the total rate of detection (typically written in units of 
events/day/kg/eV) can be calculated using the following expression 
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^ = '-^td\j{v,w4^, (28) 

where p-^ is the DM mass density and r/e is the valence electron number density per unit mass of the 
detector. For the DM density we use the value p^ = 0.3 GeV/cm'^ Q. The function /(%) is the DM 
velocity distribution in the lab frame, meaning that we need to account for the Earth's average velocity 
ve through the galaxy. We use the conventional Maxwellian distribution, truncated at a maximum escape 
velocity of v^sc (in the average rest frame of the galaxy), 

/K) = iexp(^<i±^V (29) 



where k is a normalization factor chosen such that J d^v^ /(%) = 1; 



{nv^,f/^ 



erf ( ^^ ) - A!^e-Lc/-g 



(30) 



We use the following values for the velocity parameters: average Earth velocity ve = 240 km/s, 
average DM velocity vq = 230 km/s, and DM escape velocity Vesc = 600 km/s [^. These values give us 
k = 2.504 X 10"^ (in units with c = 1). 

Performing the angular integrals (of which our cross-sections are independent), we obtain the following 
expression 



dR p^ rjf, ttUq 



^ ^ dER J^ . ^ ^ dER 



dER m^kvE 
where Vmin is the minimum DM velocity necessary to induce an interaction of recoil energy Er, 



(31) 



^2{Er + Eb) 



(32) 



4. Sensitivities 

4-1- Approach 

Using these detection rates, we can estimate the possible sensitivity for a detector such as CDMSLite. 
We assume a flat background rate of 1 event/day/kg/keV, which was provided as an experimental estima- 
tion of the expected background Q. This background rate is due to residual radioactivity, which produces 
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high energy gamma rays that can Compton scatter in the detector. To determine if the experimental reach 
is hmited by the background, we also include the sensitivity possible using germanium with no background. 
Due to the small number of events, we find that reducing the background beyond our estimate does not 
have a substantial effect on the exclusion sensitivity. 

The proposed detection method involves the measurement of single electron- hole pairs, assuming that 
the energy deposited in the initial recoiling electron will prompt the formation of secondary electron-hole 
pairs. This method of detection will limit the energy resolution to the average energy per electron-hole 



pair, which is approximately 3 eV |@, [Sj. We do not consider this process in detail, and instead focus on 
the initial deposition of energy into a single electron. The possible interactions we consider are all peaked 
very strongly at low Eji ~ eV, falling off quickly with growing recoil energy. Because of this, we focus solely 
on recoils with the lowest energy {Er < 9 eV) where the signal is most competitive with the background. 
For CDMSLite, we assume the detection setup of CDMS II p^, with 4.4 kg of germanium and 1.1 kg of 



silicon. We also consider the sensitivity of each material separately for the sake of comparison, rather than 



combine data for a total exclusion limit. Using the approach of |39|, we can find the 95% confidence level 
sensitivity possible after an experimental run of one year. 

For comparison, we also include the possible sensitivity of measuring electron recoils in a xenon-based 
detector. This approximate sensitivity is calculated using the ionization cross-sections derived earlier, 
combined with atomic data from [^^. For simplification, we only consider interactions with the valence 
electrons in the 5s- and 5p-states, which dominate the overall cross-section, and use Z^ff = 1 for the Som- 
merfeld enhancement. We assume a detector mass of 1 kg, a runtime of 1 year, and negligible background, 
in order to find the maximum possible reach of noble gas detectors. 

4-2. Results 

Results are shown in Figure Q for dipole moments and in Figure || for broken U{1) models, with 
sensitivities for silicon, germanium, and xenon detectors. As mentioned before, we consider the mass range 
600 keV < m^ < 10 GeV. The upper limit is due to current nuclear recoil experiments, most of which 
provide substantial exclusions down to ttt,-^ ~ 10 GeV Q. The lower limit is simply due to the limited 
energy available for recoils at such light masses, which suppresses the possible signal. 

For each interaction, we show the strongest current experimental and astrophysical constraints on 



the relevant parameter space. There are also detailed constraints placed by supernova cooling |41, 42 1 
and BBN [^] on DM with mass m^ < 10 MeV, but the full calculation of those constraints for these 
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Figure 4: Exclusion sensitivity at 95% confidence level possible after 1 year, for (a) electric and (b) magnetic dipole moments. 
The solid lines assume a background of 1 event/day/kg/keV, while the dashed lines assume no background. Areas above 
the curves for germanium (red), silicon (blue), and xenon (brown) would be excluded. Regions in gray are already excluded 
for all models of DM by other experiments or astrophysical data. Masses to the left of the dashed black line are potentially 
constrained by supernova cooling and BEN. While a detailed calculation of these constraints on lighter masses is beyond the 
scope of this work, it is unlikely the entire region is fully excluded. 

particular models is beyond the scope of this paper. For this reason, we simply indicate in our plots the 
mass value below which these additional bounds potentially apply. As a conceptual reference for the plots, 
the minimum dy. and fi^ values excluded by germanium without background for rriy. = 10 MeV correspond 
to {av) w 10"^^ cm^. 

As discussed in the beginning of section ^, the larger binding energy present in xenon necessitates a 
larger momentum transfer. This results in an increased form factor suppression, reducing the experimental 
reach of noble gases in comparison with that of semiconductors, as seen in both Figures |^ and ^. 

For each of these exclusion limits, germanium provides somewhat weaker limits. This is caused by our 
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Figure 5; Exclusion sensitivity at 95% confidence level possible after 1 year, for eflfective C/(l) coupling j A = ^\j -^ ) with 



(a) mA ~ 10 MeV and (b) m^ = 1 meV. The solid lines assume a background of 1 event/day/kg/keV, while the dashed lines 
assume no background. Areas above the curves for germanium (red), silicon (blue), and xenon (brown) would be excluded. 
Regions in gray are already excluded for all models of DM by other experiments or astrophysical data. Masses to the left of the 
dashed black line are potentially constrained by supernova cooling and BBN. While a detailed calculation of these constraints 
on lighter masses is beyond the scope of this work, it is unlikely the entire region is fully excluded. 



assumption of the same background per unit mass for both materials, which places a stronger restriction 
on germanium, due to its heavier nuclear mass. Germanium's reach is also slightly reduced by its smaller 
effective electron mass in the conduction band, but does become more competitive with silicon at lower 
DM masses, due to its smaller band gap. 

For EDM and MDM, the strongest general constraints on lighter masses come from colliders such 
as LEP Q]. The parameter space for larger masses (> 1 GeV) is also probed by the direct detection 
experiments XQC [||, ||], CRESST |^, CDMS [||], and XENONIO ^. Direct detection rates for 
MDM are suppressed by the extra factor of v'^ relative to EDM limits. This effect applies to both current 
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Figure 6: Detection rates (in events/eV) after 1 year, for m^ = 100 MeV. These rates assume detector masses of 4.4 kg for 
germanium (red) and 1.1 kg for silicon (blue), (a) Electric dipole with d-^ — 10~^ TeV~^. (b) Magnetic dipole with /i^ = 10^^ 
TeV"^ (c) Heavy broken U{1) with niA = 10 MeV and A = 10"''. (d) Light broken [/(I) with tua = 1 meV and A = 10~". 



nuclear recoil experiments and CDMSLite, while collider limits do not suffer from this suppression 



Current dark matter annihilation searches such as Fermi LAT [51], and bounds on dark matter anni- 



hilation rates in the early universe from WMAP and ACT |52] also provide limits to DM dipole moments. 
However, these results can only constrain symmetric DM, and do not apply to asymmetric models of DM. 
The resulting constraints are strong for MDM interactions, but the case of EDM annihilations is much 
more suppressed, making those irrelevant for our purposes. 

As discussed in section ^, models which generate dipole moments will also generate pointlike four- 
fermion interactions. Current bounds on these pointlike operators placed by collider experiments or current 
direct detection results could then potentially constrain dipole moments, as well. However, the charged 
particles which generate dipole moments will generically not couple directly to electrons or quarks. The 
generation of an effective vertex between DM and electrons would therefore require the exchange of a 
Z boson, in addition to the loop of charge intermediaries. The resulting effective operator would be 
substantially suppressed by a loop factor, the intermediary mass, and the Z boson propagator. Models 
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which generate the dipole moments considered in this paper (A > 1 TeV) would therefore be unconstrained 
by current bounds on pointhke four-fermion operators. However, the coupling of light DM to Z bosons 



which contributes to these operators is constrained by measurements of the Z width |53]. If the resulting Z 
dipole moments are approximately the same order as the EDM or MDM, we estimate the dipole moment 
scale would be constrained to be > 2 TeV. These EDM and MDM bounds are very model-dependent, and 
therefore are not included in Figure ^ 

As mentioned earlier, there are two limiting cases for a broken U{1). The first corresponds to a heavy 
mediator, for which we consider ruA = 10 MeV, and the second to a light mediator, for which we use 1 
meV. These sample masses were selected for comparison with [0. 



The strongest constraints on A, defined previously in equation (20), come from a combination of 
astrophysical data, due to the fact that A depends on both the DM self-coupling and the ^' -photon 
mixing. Details on the U{1) mixing constraints can be found in [54|, while the strongest constraints on 



DM self-coupling are explained in |g, 55, 56 1. The larger mass regions are also limited by the same direct 
detection experiments as dipole moments. 

In Figure |6|, we also provide example detection rates for each model considered above, assuming a DM 
mass of 100 MeV, as well as the same runtime (1 year) and detector masses (4.4 kg for Ge, 1.1 kg for 
Si). In the case of larger background, annual modulation in these rates due to the Earth's motion around 
the Sun would provide an important confirmation of a potential DM signal. We estimate that the annual 
modulation would be approximately 6% of the baseline average for the case of both EDM and MDM, 
which is larger than the traditional modulation associated with a velocity-independent cross-section. This 
difference arises due to the initial electron wavefunction, which provides substantial g-dependence (and 
therefore velocity-dependence) in the final cross-section. 

5. Conclusions 

We have examined the potential direct detection reach of semiconductor-based experiments such as 
CDMSLite. Experimental sensitivity to the production of single electron-hole pairs can dramatically im- 
prove the detectable energy range over that of traditional nuclear recoil methods. This enhanced energy 
range kinematically allows electron recoils to probe the parameter space of light dark matter, a region 
which remains largely inaccessible to nuclear recoils. The small energy gap present in semiconducting 
bandstructure also provides materials such as silicon and germanium with a substantial advantage rela- 
tive to noble gases. Motivated by these prospects, we have considered the possible interactions between 
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electrons and dark matter, many of which fit naturahy into weak-scale extensions to the standard model. 
We have found that semiconductor detectors are sensitive to a large range of uninvestigated parameter 
space, specifically for interactions such as dipole moments, which are enhanced at low recoil energies. Such 
dipole moments are naturally generated in many extensions of the standard model and are generically 
expected in models where the dark matter carries a conserved charge, such as asymmetric dark matter 
|T5HT9[]. Electromagnetic moments provide a unique glimpse into higher-scale physics, and we have found 
that CDMSLite can extend our current reach by orders of magnitude, up to scales as large as 10^ TeV. Light 
dark matter with dipolar or new gauge interactions remains a well-motivated alternative to the traditional 
heavy WIMP scenario, and CDMSLite would present a substantial opportunity to explore this possibility. 
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APPENDIX 

Here we present a more detailed description of ionization by dark matter. Our interaction system 
consists of the DM (x^), electron (xe), and nucleus (xn). Since the electron and nucleus originally form a 
bound state, we can change coordinate systems to more clearly show that this is a 2- to 3-body scattering 
process (bound atom + DM — )• electron + nucleus + DM). These coordinates are similar to the relative 
and center-of-mass coordinates used in strictly two-body systems. First we define new coordinates for the 
atomic electron-nucleus system, 



Xa — Xg X]\[, 

me ^ ^ ruN 

XA = ■ Xe H ■ XN- 

rUe + rriN me + m,^ 



(33) 



We can then use these to define new coordinates for the full DM-atom system. 



Xr — X^ XAi 

m,y ^ ruA 

XR = x^ H ■ XA- 

m,y. + m-A 



(34) 



These new coordinates make it simple to write the system's initial state, as well as the relative velocity 
necessary for our scattering equation, 

18 



^>,{Xa,Xr,XR) = ^{Xa)e'^^-''^e'^^-^^, 



Pr 
Vrel = 



Px 



PA 



(35) 



while the final state can easily be written in either coordinate system 



^f{xa,Xr,XR) = e^^--^-e'^^-'^'-e'^^-^^ = e'^-'^-e'^^-'^^ e'^^'^^. (36) 

Now we should be more careful and worry about normalizing these free state wavefunctions, but 
whatever normalization factors we use now will cancel out and not affect our final interaction cross-section, 
so we will ignore them. We can now use these coordinates to simplify the matrix element for the electron 
recoil, assuming the interaction Hamiltonian H is only a function of the relative displacement between the 
DM and electron 



\{f\H\t)\' = {27:Y5\kf-pi) 



'<P[ka + -^{kr -Pr 
nie 



H{k,f — Pr) 



(37) 



where H{'p) can be found by simply taking the nonrelativistic limit of the interaction amplitude Ai{p), 
with all spinors normalized to 1. The last assumptions we need to make to simplify our final cross-section 
formula is that the incoming momentum of the nucleus is zero and m-e < < mN ■ We can then integrate over 
the outgoing nuclear momentum and convert back to the more intuitive physical coordinates, resulting in 



da 



\Vrel\ 



n 



(fkf_ 

(2^ 



1 d^k^ d^ke 



27TS{Ef-E,)\{f\H\i) 



2tt S{Ef — Ei) '>p{Px — k^ — ke 



(38) 



H{P^ 



\v^\ (27r)3 (27r)3 

In the limit mjv — ^ oo, the fixed nucleus breaks translation invariance, and therefore momentum 
conservation. Perhaps a more straightforward interpretation is that the nucleus can absorb any finite 
momentum at negligible energy cost, due to the mjv suppression in its kinetic energy. Either way, energy 
conservation is now the only constraint on the electron-DM system. 

We can further simplify our final expression by taking the limit k^ « q, where q is the magnitude of 
the momentum transfer {q = Px~ k^)- We can check that this limit is valid by comparing the expressions 
for ke and the minimum possible q for a given incoming v^, 
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I (39) 

Qmin = m^v^ - \Jm\v\ - 2m^{ER + Eb)- 

We then see that q-min can be further reduced either by increasing m-^ or v^ or by decreasing Er or Er- 
If we then consider the most extreme case, with m^ = 10 GeV, v^^ = Vesc, Er = 1 eV, and Er = 1 eV, we 
find that ke ~ Qmin ~ 1 keV. However, the density of states for such a smah binding energy is neghgible, 
especially in the case of sihcon. Increasing the binding energy to values with higher detection efficiency will 
only increase the ratio of qmin to ke, as will decreasing rriy. or v^. For the bulk of the relevant parameter 
space, then, ke will be substantially smaller than even the minimum possible momentum transfer, making 
this approximation valid. 

As a pedagogical example, we now choose a simple point-vertex interaction and use the ground-state 
hydrogen wavefunction for ip. After enforcing energy conservation and integrating over the trivial angles, 
we find the expression 

da _ Wa^kemek^m^ , ^ ^ 2^-4 



dER A47r2 



X 



j dcose^{l + a\p^-k^f)-\ (40) 



The naive expectation would be for this integral to be 0{1) at low energy recoil (and therefore low 
momentum transfer), but it turns out that there is a large suppression for most of the range of integration. 
This can be understood as a phase space suppression for any angle that deviates away from forward 
scattering (or alternatively, for any momentum exchange q that deviates away from the minimum value 
qmin = Px~ ^x)- This means that our final expressions will be dominated by terms corresponding to qmin- 
In this example, we obtain 



da ^ Sakeme ^^^ ^ ^,^^^ _ ^^^,^_3 _ ^^ ^ ^^(^^ + ^^^2y.^ 



dER 3A47r2i;2 



(41) 

SakeTrUe 1 



The second apparent form of suppression comes in the factor of nie in the numerator, rather than some 
form of atom-DM reduced mass. This can be understood by considering the limiting behavior as me — > oo. 
In this limit, the cross-section rapidly grows for recoil energies Er ~ 0, and is heavily suppressed for all 
others, converging on a delta function centered at zero recoil energy (an ionized electron with no kinetic 
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energy). This limiting behavior should have been expected to appear, so it is not surprising to see the 
factor kerrie in our final answer. 

Applying this same approach to the models considered in this paper, we obtain the following full 
cross-sections 



da 



EDM 



lQa?d^ke 



Qa'^{p^ + k^f + 15a2(p^ + k^f + 11 



6(l + a2(p^ + A:^)2)3 



2 In 



Px ^ 

Px + ^> 



+ In 



ky/) 



i + g^fa 

l + a?{py. + k^ 



6aHPx - k^Y + 15a' (Px 



6(l + a2(p^-fe^)2)3 



kx? + 11 



(42) 



da 



MDM 



37rf2 



Ga^ip^ + k^)^ + 15a2(p^ + k^f + 11 



6(l + a2(p^ + A;x)2)3 



2 In 



Px 



Px + k^ 



+ In 



l + a\p^ 



kx? 



1 + a-^iPx + k^y 



Ga^iPx 



k^)^ + 15aHp^ 



kx) 



11 



+ -{l + a\p^ 



kx?)-' 



1 



6(1 + a?{p^ 



kx)')' 



(43) 



da 



hl2\^a^ke 
U(i) v\{a^m\ - 1) 



(1 + a^{Px - kx?)ijn\ + (p^ + k^f) 

{a^m\ - 1) ^^^ V((l + «'(Px + kxf){jn\ + (p^ " k^f) 



■In 



\2\2 



+ 



a^m\ - Xy - 3(a^mji - 1)(1 + a\p^ - k^Y) + 9(1 + a\p^ - k^ 

12(l + a2(p^-A:^)2)3 

{a^m\ - If - 3{a^m\ - 1)(1 + a^jp^ + k^f) + 9(1 + a^jp^ + fc^^)^)^ 

2\3 



12(l + a2(p^ + A;^)2) 



+ 



(44) 



4a2(m^ + {py. - k^Y) Aa'^{m\ + {p^ + k^Y) _ 
If we then take the appropriate limits of these cross-sections, we obtain the approximate results given 
earlier. 
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